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Abstract 

The effects of strong magnetic fields on neutron star matter are investigated in the quark-meson 
coupling (QMC) model. The QMC model describes a nuclear many-body system as nonoverlapping 
MIT bags in which quarks interact through self-consistent exchange of scalar and vector mesons in 
the mean-field approximation. The results of the QMC model are compared with those obtained 
in a relativistic mean-field (RMF) model. It is found that quantitative differences exist between 
the QMC and RMF models, while qualitative trends of the magnetic field effects on the equation 
of state and composition of neutron star matter are very similar. 
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I. INTRODUCTION 



Recent observations of soft gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs) have suggested that the surface magnetic field of young neutron stars could be 
of order 10^^ — 10^^ G ill. On the other hand, it is estimated that the interior field in neu- 
tron stars may be as large as 10^*^ G . Motivated by the existence of strong magnetic 
fields in neutron stars, theoretical studies on the effect of extremely large fields on dense 
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matter and neutron stars have been carried out by many authors 
The inclusion of hyperons and boson condensation has also been investigated 
It was found that magnetic fields could change the composition of matter dramatically, and 
the softening of the equation of state (EOS) caused by Landau quantization is overwhelmed 
by the stiffening due to the incorporation of the nucleon anomalous magnetic moments for 
field strengths B > 10^ {B^ = 4.414 x 10^^ G is the electron critical field). So far, most 
of these studies have employed the relativistic mean-field (RMF) theory, which is a field 
theoretical approach at hadron level and has been widely used in the description of nuclear 
matter and finite nuclei 141. It is well known that various RMF models could provide similar 

n 

nuclear saturation properties but different behaviors at high density. In Ref. p[, the model 
dependence of strong magnetic field effects within the RMF theory has been investigated. 
It is instructive to study the properties of neutron star matter in the presence of strong 
magnetic fields by employing various approaches. 

In this paper, we adopt the quark-meson coupling (QMC) model for the investigation 
of strong magnetic fie d effects on neutron star matter. The QMC model was originally 
proposed by Guichon in which the quark degrees of freedom are explicitly taken into 
account. The QMC model describes a nuclear many-body system as nonoverlapping MIT 
bags in which quarks interact through self-consistent exchange of scalar and vector mesons 
in the mean-field approximation. The mesons couple not to point-like nucleons but directly 
to confined quarks. In contrast to the RMF approach, the quark structure of nucleons plays 
a crucial role in the QMC model, and the basic coupling constants are defined at quark 
level. The QMC model has been subsequently extended and a ppli ed to various problems 
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Furthermore, the 



of nuclear matter and finite nuclei with reasonable success 
model has also been used to investigate the properties of neutron stars with the inclusion of 
Kype.,., ,ua..U. and Uao. co,*,.a«o„ iJflQ. It . .te.e.,,, to extend the QMC 
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model for the study of a nuclear many-body system in strong magnetic fields, and make a 
systematic comparison between the QMC and RMF models. 

In general, a nuclear many-body problem has to be worked out in the QMC model by 
two steps. First, one can calculate the nucleon properties in medium by using the MIT bag 
model with an external scalar field a generated by nuclear medium. Then in the second 
step, the entire nuclear system is solved at hadron level with the effective nucleon mass 
as a function of a obtained in the first step. The main difference between the QMC and 
RMF models is reflected in the dependence of effective nucleon mass on scalar field a. In 
the QMC model, the scalar field plays a vital role in determining the nucleon properties 
such as effective mass and radius, while the vector fields do not cause any change of nucleon 
properties but appear merely as energy shifts. When we treat a nuclear many-body system 
in strong magnetic fields, we assume that the magnetic fields do not change the internal 
structure of nucleons, then nucleons in meson and magnetic fields could be treated quite 
similarly as in the RMF models. With the effective nucleon mass obtained at quark level, 
we can investigate the properties of neutron star matter in strong magnetic fields. 

This paper is arranged as follows. In Sec. II, we briefly describe the QMC model for 
neutron star matter and present relevant equations in strong magnetic fields. In Sec. Ill, we 
show and discuss the numerical results in the QMC model and make a systematic comparison 
with the results of the RMF model. The conclusion is presented in Sec. IV. 

II. QUARK-MESON COUPLING MODEL FOR NEUTRON STAR MATTER IN 
STRONG MAGNETIC FIELDS 

The QMC model describes a nuclear many-body system as nonoverlapping spherical bags 
in which quarks interact through self-consistent exchange of scalar and vector mesons in the 
mean-field approximation. To perform the calculation for neutron star matter in strong 
magnetic fields, we first study the nucleon properties with external meson and magnetic 
fields by using the MIT bag model. The fields in the bag are in principle functions of 
position, which may cause a deformation of the nucleon bag. For simplicity, we neglect the 
spatial variation of the fields over the small bag volume, and take the values at the center 
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of the bag as their average quantities Then the quarks in the bag satisfy the Dirac 
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equation 

' '(1 + 3-3). .,1^,.0. (1) 



where g^, g^, and g'^ are the quark-meson couphng constants, and rriq is the current quark 
mass, cr, ujfj^, p3^, and are the values of the meson and magnetic fields at the center of 
the nucleon bag. The normalized ground state for a quark in the bag is given by 

V,,(r.<)=^,e-.'/«( (2) 

ytf3ga-rji{xgr/R) J V47r 

where 



.Vtq — Rm* 



Af-' = 2R'f,{xg) - 1) + Rml/2] /x] , (4) 



with Vtq= -y/a^q + {Rm*y, m* = niq + g^a. R is the bag radius, and Xq is the quark spinor. 
The boundary condition, io{xq) = Pqji{xq), at the bag surface determines the eigenvalue 
Xq. The energy of a static nucleon bag consisting of three ground state quarks is then given 
by 

^bag = 3^ - I + l^R'B,ag , (5) 

where the parameter Z accounts for zero-point motion, and -Bbag is the bag constant. The 
effective nucleon mass is then taken to be 

= Ebag • (6) 

The bag radius R is determined by the equilibrium condition dM'^/dR = 0. In the present 
calculation, we take the current quark mass = 5.5 MeV. The parameter = 210.854 
MeV and Z = 4.00506 are determined by reproducing the nucleon mass in free space M^v = 



939 MeV and the bag radius R = 0.6 fm as given in Ref. P0|. We note that uj and p 
mean fields appear merely as the energy shifts which do not cause any change in nucleon 
properties, therefore the effective nucleon mass obtained in the QMC model depends on the 
a mean field only. 

For neutron star matter consisting of a neutral mixture of neutrons, protons, electrons, 
and muons in f3 equilibrium in the presence of strong magnetic fields, the total Lagrangian 
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density at hadron level can be written as 



b=n,p 

1„ 1 
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l=e,fi 
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where ipb and ipi are the nucleon and lepton fields, respectively. = (0, 0, Bx, 0) refers 
to a constant external magnetic field B along the 2;-axis. The field tensors for the u, p, and 
magnetic field are given by = d^ujy-dyU^, Ri^„ = d^pi^-dypi^, and F^y = d^A^-dyA^. 
It has been argued in Ref. jj] that the contributions from anomalous magnetic moments 
of nucleons should be considered for large field strengths. Here the anomalous magnetic 
moments of nucleons are included with Hp = p^ i.Qp/'^ ~ 1) = 1.7928 /iTv and /t„ = p^Qn/^ = 
— 1.9130 pn, where pn is the nuclear magneton. As for leptons, the effects of anomalous 
magnetic moments of electrons and muons on the EOS are very small as shown in Ref. j^, 
and the electron anomalous magnetic moments could be efficiently reduced by high-order 
contributions from the vacuum polarization in strong magnetic fields j^]. Therefore, we 
neglect the anomalous magnetic moments of leptons in the present work. The coupling 
constants at hadron level, and gp, are related to the corresponding quark-meson couplings 
as gi^ = 3g^ and gp = g'j^ j2^. The quark-meson coupling constants gl = 5.95 7, q'^. = 2.994, 
and = 4.325 are determined by fitting the properties of nuclear matter j20], and the 
meson masses mo- = 550 MeV, = 783 MeV, and irip = 770 MeV are used in the present 
calculation. 

From the Lagrangian density given in Eq. (|7j), we obtain the following meson field equa- 
tions in the mean-field approximation 

mla = -^ip^,+p:), (8) 
mluJo = gu.{pP, + p",) , (9) 
"ijpso = 9p{ff, - p") , (10) 

while the Dirac equations for nucleons and leptons are given by 

{ilpd^^ - qbl^,A^' - M*^ - g^'y^o - Qpl'^r^Pm - l^^b^x^.i^'^") V-fe = , (11) 
{ili.d^ - qa^.A^' -mi)tPi = 0. (12) 



The energy spectra for protons, neutrons, and leptons (electrons and muons) are given by 

'kl+[ \/M*^^+2uqpB - sKpB ] + g^uo + gpPso , (13) 

(14) 
(15) 



E: = Jk^,+ {^^M*^ + kl + kl - sKnB^ + g^uo - g^pm , 



E. 



'kl + mf+2u\qi\B 



where v = n + 1/2 — sgn (q) s/2 = 0, 1, 2, . . . enumerates the Landau levels of the fermion 
with electric charge q. The quantum number s = ibl are for spin-up and spin-down cases. 
As discussed in Ref. jsf, the Dirac spinors are no longer eigenfunctions of the spin operator 
along the magnetic field direction when the anomalous magnetic moments are taken into 
account. However, s reduce to the eigenvalue of cr^ as the anomalous magnetic moment tends 
toward zero. We note that the first Landau level to be occupied is the state with z/ = 0, 
s = -|-1 for the proton and u = 0, s = —1 for the lepton, while the neutron first occupies 
the state with s = — 1 due to its negative anomalous magnetic moment. This argument is 
consistent with the results in nonrelativistic approach P|. 

The expressions of the scalar and vector densities for protons and neutrons are given by 



qpBM*^ ^ ^ y/Mj^^ + 2uqpB - sk^B 



2tx' 

V s 

27,^ 5Z 5Z ^L'^,s ' 



^/M*fT2^ 



^/M^fT2uq^ - SKpB 



47r2 ^ 



kl^E'f-iM^^-SK^Bfln 



M 



N 



SK„B 



(16) 
(17) 
(18) 



- SKnB TT 



+Ef ^arcsm 



(19) 



where k^^^^ and k^ ^ are the Fermi momenta of protons and neutrons, which are related to 
the Fermi energies E^ and E^ as 

'M*^^ + 2uqpB - SKpB^ , (20) 



Ef 
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ky, + (M^ - sKnBy 



(21) 
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The chemical potentials of nucleons and leptons are given by 



(22) 
(23) 
(24) 



For neutron star matter with uniform distributions, the composition of matter is deter- 
mined by the requirements of charge neutrality and /3-equilibrium conditions. In the present 
calculation, the /5-equilibrium conditions are expressed by 



(25) 
(26) 



and the charge neutrality condition is given by 



p'v = P'v + Pv 



(27) 



where the vector density of leptons has a similar expression to that of protons 



(28) 



We solve the coupled Eqs. (injl- lfTUj) . (P^. (pUj). and (P7j) self-consistently at a given baryon 
density = p^ + p" in the presence of strong magnetic fields. The energy density of neutron 
star matter is given by 



111 

= £p + + + + 2"^^^^ + -j^l^^l + 2"^P^30 ' 

where the energy densities of nucleons and leptons have the following forms: 



(29) 



47^2 



EE 



+ ( \I^N + 2i^gp5 - sKpB ) In 



E {^^/.^f - y-nBEf (^arcsin 



^M*^ + Ivq^B - SKpB 



,(30) 



— SKnB TT 



SKnB — SKnB 



X 



£1 



{M*^ - SKnB) k%E] + {Ml, - SKnBfln 

qi\B 



47r2 



EE 



k^f,^,sEf+ {m^+2u\qi\B) In 



M*^ - SKnB 



^ymJ+2u\qi\B 



(31) 
(32) 
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The pressure of the system can be obtained by 

Pm = f^iPv ~ = f^npB ~ £m, (33) 

i 

where the charge neutrahty and /?-equihbrium conditions are used to get the last equahty. 
Note that the contribution from electromagnetic fields to the energy density and pressure, 
Ef = Pf = -B^/Stt, should be taken into account in the calculation of the EOS. 



III. RESULTS AND DISCUSSION 



In this section, we present the properties of neutron star matter in strong magnetic fields 
by using the QMC model. The effective nucleon masses in the QMC model are obtained 
self-consistently at quark level, which is the main difference from the RMF model. The 
comparison of results in several RMF models has been investigated extensively in Ref. P|. 
It was found that quantitative differences persist between those models, while qualitative 
trends of magnetic field effects are very similar. Here we would like to make a systematic 
comparison between the QMC and RMF models, which are based on different degrees of 
freedom. 

In Fig. we show the matter pressure Pm as a function of the matter energy density Em 
for the field strengths B* = B/B^ = 0, 10^ and 10^ The resuhs of the QMC model in the 
upper panels are compared with those of the RMF model with TMl parameter set in 
the lower panels. In Ref. P], the RMF model with GM3 j^l and several other parameter 
sets has been used and compared with each other. Here we prefer to make a comparison 
with the TMl model, which includes nonlinear terms both for a and u mesons. The TMl 
model was determined in Ref. 2J| by reproducing the properties of nuclear matter and finite 



nuclei inclu ding neutron-rich nuclei, and has been widely used in many studies of nuclear 
physics We show the results with and without the inclusion of nucleon 



anomalous magnetic moments in the right and left panels, respectively. In both QMC 
and TMl models, the Landau quantization of charged particles causes a softening in the 
EOS as shown in the left panels, while the inclusion of anomalous magnetic moments leads 
to a stiffening of the EOS as shown in the right panels. These effects become significant 
as magnetic fields increase above B* ~ 10^. The softening of the EOS due to Landau 
quantization can be overwhelmed by the stiffening due to the incorporation of anomalous 
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magnetic moments with increasing magnetic fields for B* > 10^. Tlie EOS in tlie QMC 
model is sliglitly stiffer tlian tfie one in tlie TMl model. We show in Fig. |21 the effective 
nucleon masses as a function of baryon density ps again for B* = 0, 10^, and 10^. It is found 
that the influence of magnetic fields on the effective masses is not observable until B* > 10^. 
Without the inclusion of anomalous magnetic moments, the effective masses for B* > 10^ 
are smaller than the field-free values, but they become to be larger than the field-free values 
in strong magnetic fields when the anomalous magnetic moments are included. It is clear 
that the changes of effective masses in the QMC model are much smaller than those in the 
TMl model. 

It is instructive to discuss the spin polarization of nucleons and electrons in neutron star 
matter with the presence of strong magnetic fields. For charged particles such as electrons 
and protons, the energy spectrum is characterized by the quantum number of Landau level 
u and the spin projection s. If we ignore the anomalous magnetic moments of the charged 
particles, all energy levels except the ground state are doubly degenerate with opposite spin 
projections. The ground state of electrons (protons) is a single state with u = 0, s = —1 
{s = +1). As the field strength increases, the energy gap between Landau levels becomes 
larger and larger, and at a critical field, all particles only occupy the first Landau level 
with u = 0, s = —1 for electrons and s = +1 for protons. This is called complete spin 
polarization which occurs as the energy gap between Landau levels is comparable to the 
Fermi momentum of the particle. It is the order of B* ~ 10^ for electrons and protons 
in neutron star matter at the densities of interest. Note that neutrons can not be altered 
by external magnetic fields when the anomalous magnetic moments are ignored. It is very 
interesting to estimate the contributions from anomalous magnetic moments on the spin 
polarization of nucleons. When the anomalous magnetic moments are included, the double 
degeneracy of Landau levels for protons disappears, and the energy spectra for both neutrons 
and protons are spin-dependent. The ground state of protons (neutrons) is spin- up (spin- 
down). The next level to be occupied for protons is spin-up {u = l,s = +1), then the 
third level may be spin-down (z/ = 1, s = —1) or spin-up {u = 2, s = +1), which should be 
determined by comparing their energy spectra. Neutrons can be completely polarized with 
the inclusion of anomalous magnetic moments. The spin-up states begin to be occupied by 
neutrons for ^ 4 M'^B, so complete spin polarization of neutrons appears at the 
order of B* ~ 10^. We define the proton spin polarization as the ratio of spin-up proton 
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density to total proton density, p^(s = The results with and without the inclusion 

of anomalous magnetic moments are shown in the middle and left panels of Fig. El while 
the neutron spin polarization defined as the ratio of spin-down neutron density to total 
neutron density, p"(s = — l)/p", is shown in the right panels. It is obvious that both proton 
and neutron spin polarizations increase with decreasing baryon density and increasing field 
strength, and a complete spin polarization occurs at a critical B* for a given baryon density. 
By comparing the left and middle panels, one can find that the inclusion of anomalous 
magnetic moments leads to a significant increase in the proton spin polarization. We note 
that there is no noticeable difference between the QMC and TMl models in Fig. IHl 

In Fig. m the proton fraction Yp = p^/ps and the muon fraction = p'^/ps are plotted 
as functions of the field strength B* for the baryon densities pb = 0.075, 0.15, 0.30, and 
0.60 fm~^. One can easily obtain the neutron fraction Yn = P^/pb and the electron fraction 
Ye = Pvl Pb by the relations: y„ = 1 — 1^ and Ye = Yp — Y^. The results with and without 
the inclusion of anomalous magnetic moments are shown in the right and left panels. It is 
obvious that the composition of neutron star matter depends on both the baryon density 
Pb and the magnetic field strength B* . When the external magnetic fields are absent, the 
composition of matter is dominated by neutrons with a small proton fraction. It is clear that 
the proton fraction remains unaltered from the field-free case for relatively low fields. As the 
field strength increases to the value of having complete proton spin polarization, a significant 
enhancement of the proton fraction is observed due to large reductions in chemical potentials 
of protons and leptons caused by strong magnetic fields, and then protons become to be 
dominant for B* > 10^. On the other hand, lepton fractions increase with the increasing of 
proton fraction by the requirement of charge neutrality. The effect of strong magnetic fields 
leads to a reduction in the electron chemical potential, so that muons disappear when the 
electron chemical potential is less than the rest mass of the muon. It is evident that the 
inclusion of anomalous magnetic moments causes a slight increase in the proton and lepton 
fractions. 

IV. SUMMARY 

In this paper we have investigated the effects of strong magnetic fields on the properties 
of neutron star matter within the QMC model. The QMC model describes a nuclear many- 
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body system as nonoverlapping MIT bags in which quarks interact through self-consistent 
exchange of scalar and vector mesons in the mean-field approximation. The effective nucleon 
masses in the QMC model are obtained self-consistently at quark level, which is an important 
difference from the RMF models. It is clear that the effects of strong magnetic fields become 
significant only for field strength B* > 10^. The Landau quantization of charged particles 
causes a considerable enhancement of the proton fraction and a softening in the EOS, while 
the inclusion of nucleon anomalous magnetic moments leads to a stiffening of the EOS. With 
the inclusion of anomalous magnetic moments, a complete spin polarization of fermions 
occurs at a critical B* which depends on the baryon density. The complete spin polarization 
of neutrons appears at higher field than the one of protons. 

Wc have made a systematic comparison between the results of the QMC and RMF models. 
It is found that quantitative differences exist between these models, while qualitative trends 
of magnetic field effects on the EOS and composition of neutron star matter are very similar. 
In the interior of neutron stars hadronic matter may undergo a phase transition to deconfined 
quark matter. To estimate the phase transition in the presence of strong magnetic fields, 
the QMC model should be applicable since it describes hadronic matter in terms of quark 
degrees of freedom. 
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FIG. 1: (Color online) The matter pressure Pm versus the matter energy density Em for the field 
strengths B* = 0, 10^, and 10^. The results of the QMC model with and without the inclusion 
of anomalous magnetic moments are shown in the upper-right and upper-left panels, respectively. 
Those in the TMl model are also shown in the lower panels for comparison. 
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FIG. 2: (Color online) The effective nucleon mass M]l^/Mjv as a function of the baryon density 
Pb for the field strengths B* = 0, 10^, and 10®. The results with and without the inclusion of 
anomalous magnetic moments are shown in the right and left panels, respectively. 
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FIG. 3: (Color online) The spin polarization of protons (left and middle panels) and neutrons (right 
panels) as functions of magnetic field strength. The dotted (blue), solid (black), dashed (green), 
and dashed-dotted (red) lines correspond to pB = 0.075, 0.15, 0.30, and 0.60 fm~^, respectively. 
The neutron spin polarization is purely due to the interaction of anomalous magnetic moments 
with the magnetic field. The inclusion of anomalous magnetic moments leads to an increase in the 
proton spin polarization. 
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FIG. 4: (Color online) The proton and muon fractions as functions of magnetic field strength. 
The dotted (blue), solid (black), dashed (green), and dashed-dotted (red) lines correspond to 
Pb = 0.075, 0.15, 0.30, and 0.60 fm^^, respectively. The neutron and electron fractions could be 
obtained by y„ = 1 — 1^ and Ye = Yp — Y^. 
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